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ABSTRACT 



5^ ' Tidal interaction in a coalescing neutron star binary can resonantly excite the g-mode oscillations of 

, the neutron star when the frequency of the tidal driving force equals the intrinsic g-mode frequencies. We 

iy-j , study the g-mode oscillations of cold neutron stars using recent microscopic nuclear equations of state, where 

we determine self-consistently the sound speed and Brunt- Vaisala frequency in the nuclear liquid core. The 
properties of the g-modes associated with the stable stratification of the core depend sensitively on the 
pressure-density relation as well as the symmetry energy of the dense nuclear matter. The frequencies of the 
first ten g-modes lie approximately in the range of 10 — 100 Hz. Resonant excitations of these g-modes during 

o 

lag develops even in the absence of fluid viscosity. However, since the coupling between the g-mode and the 
tidal potential is rather weak, the amount of energy transfer during a resonance and the induced orbital 



the last few minutes of the binary coalescence result in energy transfer and angular momentum transfer from 
the binary orbit to the neutron star. The angular momentum transfer is possible because a dynamical tidal 



2 ' phase error are very small. 

Resonant excitations of the g-modes play an important role in tidal heating of binary neutron stars. 
Without the resonances, viscous dissipation is effective only when the stars are close to contact. The resonant 
oscillations result in dissipation at much larger orbital separation. The actual amount of tidal heating depends 
on the viscosity of the neutron star. Using the microscopic viscosity, we find that the binary neutron stars 
are heated to a temperature ~ 10 s K before they come into contact 
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1. INTRODUCTION 

Coalescing neutron star-neutron star (NS-NS) and neutron star-black hole (NS-BH) binaries are believed 
to be one of the most promising sources of gravitational radiation (e.g., Schutz 1986, Thornc 1987) that 
could be detected by the new generation of instruments such as the Laser Interferometer Gravitational- Wave 
Observatory (LIGO, see Abramovici et al 1992) and its French-Italian counterpart VIRGO (Bradaschia et 
al 1990). Extrapolation of the estimated galactic formation rates of compact binaries based on the observed 
local population of binary radio pulsars leads to an estimate of the rate of coalescence in the Universe of 
about 1CP 7 yr _1 Mpc -3 (Narayan, Piran & shemi 1991, Phinney 1991). Extracting gravity wave signals 
from noise requires accurate theoretical waveforms in the frequency range 10 — 1000 Hz, corresponding to 
the last few minutes of the binaries' life (e.g, Cutler et al. 1993). Effects of tidal interaction on the orbital 
decay and gravitational radiation waveforms have been found to be small except during the final stage of 
coalescence (Kochanek 1992, Bildsten & Cutler 1992, Lai, Rasio & Shapiro 1993, 1994a, 1994b). However, 
all these studies implicitly assume equilibrium or quasi-equilibrium tides, while dynamical tidal responses of 
the stars are completely neglected. Such dynamical tides are possible when the interaction potential between 
the two components of the binary resonantly excites the intrinsic oscillations of the neutron stars. 

Recent studies of neutron star oscillations have demonstrated the existence of low frequency g-modes 
of various kinds. McDermott, Van Horn & Scholl (1983) and McDcrmott, Van Horn & Hansen (1988) 
calculated the g-modes associated with thermally induced buoyancy for warm (hot) neutron stars. Finn 
(1987) considered crustal g-modes due to the composition discontinuities in the outer envelope of cold neutron 
stars. Reisenegger & Goldreich (1992, hereafter RG1) studied g-modes associated with the buoyancy induced 
by proton-neutron composition gradient in the cores of neutron stars. Regardless of the nature of the g-modc 
oscillations and various theoretical uncertainties, the existence of such g-modes indicates that a neutron star 
in a binary system can be resonantly excited when the frequency of the tidal driving force equals the g-mode 
frequencies. These resonant oscillations and their consequences to the binary orbital decay and the neutron 
stars themselves are the subjects of this paper. 

For a star in a binary system, the tidal potential due to the companion acts like an external perturbing 
force, with a driving frequency 2AS1 = 2(£1 — fl s ), where fi and il s are the orbital angular frequency and the 
spin angular frequency of the star (the factor of 2 results from the quadrupole nature of the tidal potential, 
and will become apparent in §2). The star behaves as a collection of harmonic oscillators with different eigen- 
frequencies. Let us consider just one such oscillation mode with frequency w. Schematically, the equation 
governing the tidal distortion £ can be written as 

£ + w 2 £- — C°ce 2iAot , (1.1) 

tyisc 

where t visc is the timescale for viscous dissipation (or other dissipations) in the star. For stationary state, 
we have 

? lo 2 - 4Af2 2 - 2iAn/t msc ' { ' 
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Since in general the viscous dissipation timescale is much larger than the dynamical timescale of the star, 
we see that when 2A0 = u, the amplitude of £ can be very large. 

Such resonant oscillations in close binaries were first considered by Cowling (1941) in his seminal paper 
on nonradial stellar oscillations. He noted that the coupling between the tidal potential and the high-order 
g-mode is weak, therefore a close resonance is needed to greatly enhance the equilibrium tide (see also Zahn 
1970). He also suggested that such close resonance will be destroyed by second order effects due to the large 
horizontal displacement of the high-order g-mode. Clearly, for long-lived binary systems, such as binaries 
consisting of main sequence stars, what exactly happens to the stars at resonance is a complicated issue 
(e.g., Zahn 1977, Nicholson 1978, Goldreich & Nicholson 1989). 

For an inspiraling neutron star binary, the problem becomes "cleaner". The rapid decay of the orbit 
induces a dynamical tidal lag, and implies that the orbital frequency does not stay at a particular resonant 
value very long. Thus the rapid decay of the orbit "kills" the infinity inherent in equation (1.2) even without 
dissipative or other higher-order effects. Moreover, we shall see that the resonant tidal distortion is small 
and a linear approximation for the excited oscillations is adequate. 

In this paper, we shall concentrate on the excitations of core g-modes. These modes involve the bulk 
region of a neutron star, while the g-modes driven by entropy gradient and density discontinuities are confined 
to its outer envelope (Finn 1987, McDermott, Van Horn & Hansen 1988). Previous study (RG1) has yielded 
a qualitative understanding of the basic properties of the core g-modes. However, the detailed properties of 
these g-modes, especially the strengths of their coupling to the tidal field, which directly affect the amplitudes 
of the excitations, depend sensitively on the values of the Brunt- Vaisala frequency in different regions of the 
star. The calculations by RG1 were based on an approximate, inconsistent ansatz for the Brunt- Vaisala 
frequency in the core. We examine some recent microscopic equations of state of nuclear matter, and use 
them to obtain the self-consistent Brunt- Vaisala frequency, upon which our calculations of g-modes are based 
(see §4 and §5). 

A related question of interest concerns how much energy is dissipated into heat in coalescing binary 
neutron stars. Meszaros & Rees (1992) considered the consequence of such tidal heating in the context of 
cosmological Gamma ray burst models. Significant tidal heating of the neutron stars can induce mass loss 
from the stars even before merging takes place. Such mass loss can lead to baryon contamination of the 
possible Gamma ray burst emission. Clearly, this result depends on the viscosity of the neutron star matter. 
We show that the resonant excitations of g-modes are important in the viscous heating of binary neutron 
stars (§8). 

Some aspects of the problem of resonant excitations of g-modes have also recently been considered 
by Reisenegger & Goldreich (1994, hereafter RG2). They calculated the total energy transfer from the 
binary orbit to the star in a resonance using the formalism of Press & Teukolsky (1977, hereafter PT), 
who considered a similar problem in the context of parabolic encounter of two stars. We present here a 
different, more complete analysis of oscillations in coalescing binary neutron stars, including dynamical tidal 
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lag, energy transfer and angular momentum transfer due to resonances, etc. Our method enables us to study 
the time-dependence of the resonant energy transfer process (rather than the total energy transferred), as 
well as the non-resonant, quasi-equilibrium tidal effects. We include comparisons of our results with those 
of RG2. 

In §2 we develop the formalism to study tidal excitations of normal modes in inspiraling binary stars, 
where we show that the problem can naturally be divided into two parts. The dynamical part is discussed 
in §3, while the physics of core g-modes and the calculations of the mode properties are discussed in §4 and 
§5. In §6 energy transfer and angular momentum transfer from the binary orbit to the star are considered. 
The effects of such transfer on the orbital decay and the induced orbital phase error are discussed in §7. 
Tidal heating of binary neutron stars are considered in §8. We present our conclusions in §9, where we also 
discuss the importance of g-mode resonances in coalescing white dwarf binaries. 

2. DYNAMICAL TIDAL INTERACTION IN COALESCING BINARIES: 

THE FORMALISM 

Consider a star with mass M and radius R in circular orbit with its companion, which we treat as a point 
object with mass M' . When M' also has a finite size, there is additional quadrupole-quadrupolc coupling 
between M and M'; but this is a higher-order effect which will be neglected. Also, general relativistic effects 
introduce a small correction to our analysis; these effects are neglected. 

The interaction potential due to M' is 

where r is the position vector of a fluid element in star M, and D =T)(t) specifies the position of the point 
mass. In a comoving coordinate system (r, 9, <f>) with its origin at the center of M and the orbital plane at 
9 = 7r/2, we have D(£) = (D(t), ir/2, 3>(i)), and the potential can be expanded into spherical harmonics 

4-7T r l 



lm 

= -GM'^Wirn^y^e-^YUe,^), 



(2.2) 



lm 

(using the notation of PT), where the numerical coefficient Wi m is 



W lm = (-) 



(l+m)/2 



rn 



(2.3) 



Here the symbol (— ) k is to be interpreted as zero when k is not an integer. In equation (2.2), the I = and 
I = 1 terms can be be dropped, since they are not relevant for tidal deformation. Although our equations 
are valid for general I, throughout this paper we shall only consider the leading quadrupole term, I = 2, for 
which the relevant coefficients Wi m are 

1/2 



/7T\l/2 1 /67r\ J 

W 20 = -(-) , W 2 ±1=0, T^ 2±2 = -fyj 



(2.4) 



In the linear approximation, the effect of the tidal potential on star M is specified by the Lagrangian 
displacement £ (r, t) of a fluid element from its unperturbed position. The equation of motion can be written 

as 
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P§{2+£)Z=-PVU, (2.5) 

where p is the density, and C is an operator which specifies the internal restoring forces of the star. We now 
analyze £(r, t) in to the normal modes of the star 

f(r,t)=5> a (t)£,(r), (2.6) 

a 

where a specifies the "quantum number" of a normal mode. The mode cigenfunction £ a satisfies 

(C-pu;l)£ a (r)=0, (2.7) 

where u a is the angular frequency of the eigenmode. For a spherical star, the normal modes are labeled 
by the spherical harmonic indices, I, m, and by a "radial quantum number" n, i.e., {a} — {n,l,m}. For 
convenience, the indices a and {n, I, m} will be used interchangeably throughout the paper. For a particular 
spheroidal mode, the eigenfunction can be written as the sum of the radial and tangential pieces 

&(r) = inlm{v) = [C' nl (r)e r + r^ nl (r)V]Y lm (6,<j>), (2.8) 

where e r is the unit radial vector. Note that toroidal modes, which are characterized by identically vanishing 
Eulerian variation of density (e.g., Cox 1980, p222), are not excited by potential forces such as the tidal 
interaction considered here. Using the orthogonal relations for £ Q and the normalization 



<Pxp&-Sa>=S aa ., (2.9) 



equation (2.5) can be reduced to a set of equations for a a (i): 



a a+ .la a = GM '^ Qnl e-^\ (2.10) 



where we have defined a tidal coupling coefficient 

Q nl = Jd 3 xpC lm -V[r l Y lm (0^)] 

= [ R plr l+1 dr[C nl (r) + (l + l)^ nl (r)], 
Jo 



(2.H) 



(see Zahn 1970, PT). Since £ is real, we have a* m=2 — a m =-2- 

Equation (2.10) then essentially determines the response of the star to the tidal potential, provided we 
know the functions D(t) and $(£) which specify the decaying orbit. Using the quadrupole formula for the 
gravitational radiation of two point masses, the rate of change of the orbital separation is given by 

• 64G* 3 MM'(M + M') 

5c 5 D 3 1 ' ' 
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The use of this point-mass formula is adequate, since the effect of tidal interaction on the orbit is small 
(see §6; a fully self-consistent formalism, incorporating the feedback of the tidal interaction on the orbital 
trajectory, is discussed in Lai 1994). It is convenient to define an orbital decay time to as 

D 5c 5 D± 
D \D\ 64G 3 MM'(M + M')' V ' ' 

The orbital phase function <E>(i) is 

$(£) = J dtMl = J dt(fl-n s ), (2.14) 
where f2 s is the spin of the star, and is the orbital angular frequency, for which we have 

When Q s ^ 0, the normal modes of the star become much more complicated, especially when fi s becomes 
comparable to the mode frequencies. Throughout our paper, we shall assume f2 >> fl s , so that the eigen- 
modes can be adequately approximated by those of a non-rotating spherical star, for which u a = uo n i does 
not depend on m. Kochanek (1992), and Bildstcn and Cutler (1992) have shown that the viscosity is too 
small for significant spin-up of the neutron star during the inspiral. We point out that even in the absence 
of fluid viscosity, angular momentum can be transferred to the star due to dynamical tidal lag. In §6, we 
will show that such dynamical spin-up of the neutron star is also negligible. 
It is convenient to define a function b a (t) such that 

a a (t) = {GM'W lm Q nl )b a {t)e- m ^\ (2.16) 

Equation (2.10) then becomes 

b a — 2im£lb a + (w^ — m 2 2 — im(l)b a = ■ (2-17) 

Note that apart from the mode frequency, b a (t) does not depends on the properties of the mode and the 
structure of the star. The mode properties enter through the tidal coupling coefficient Q n i. Thus the whole 
problem is divided into two parts: Find b a (t) by solving equation (2.17) for a given u a , and calculate the 
frequencies and tidal coupling coefficients Q n i for different modes. 

The natural units of mass, length, time, and energy in this problem are M, R, (i? 3 /GM) 1 / 2 , and 
GM 2 /R respectively. We shall use these units as well as real physical units throughout the paper, whichever 
is more convenient. Using the natural units to nondimensionalize Q n i, £ Q and b a , we have 

|(r,t) = h Y, WimQm { ^y +1 e- im *W£a(r), (2.18) 

(only I = 2 terms are included in the sum), where h a is simply the typical equilibrium tide height 



The equilibrium tide corresponds mainly to the f-mode distortion, for which b a ~ (R/D) 3 and Q n i ~ 1, as 
we shall see later. 



3. DIMENSIONLESS DYNAMICAL TIDAL AMPLITUDE 



In this section, we consider the dimensionless tidal amplitude b a {t) for a given mode frequency u> a - This 
determines the dynamical aspects of the tidal excitation. 

3.1 Solution Prior to The Resonance 

From equation (2.4), we see that b{t) is nonzero only for to = or m = ±2. For m = 0, the tidal 
potential is quasi-static, and equation (2.17) becomes 

1 



b a + Lo 2 a b a 



D i+i- 



(to = 0). 



(3.1) 



When the orbital decay timescale to is much longer than the oscillation period l/u> a , the function b a (t) is 
approximately equal to the quasi-equilibrium value b a . eq : 

1 



b a (t) ~ b a , eq (t) 



(to = 0). 



(3.2) 



wiD(t) ,+1 ' 

Now consider the dynamical tide (m = ±2). Resonance occurs when uo a = 2Aft ~ 2ft, at a resonance 
radius given by 



Da 

R 



4(1 + q)GM 



"I 1/3 



;r 3 



(3.3) 



where q = M'/M is the mass ratio. Before reaching this resonance, when uj 2 — (2ft) 2 » an approximate 
solution for b a (t) can be found. Dropping the b a and b a terms in equation (2.17), we obtain a zero-th order 
solution 



b^ = 



1 



D'+!(w2 - TO 2 ft 2 )' 



Thus the zero-th order term in b a is 

b n 



„.£) 2TO 2 ftf2 



£) CJ 2 -TO 2 ft 2 



b a . 



(3.4) 



(3.5) 



Substituting (3.5) into equation (2.17), and dropping the term b ~ we obtain an approximate solution 

for b a (t), valid for u? a - 4ft 2 >> |f2|: 

-l 



b a (t) ~ 



1 



£)i+l 



2 2^2 f„/, ,n A 4m 3 ft 2 ft 

UJ 2 a - TO 2 ft 2 + t [2(1 + 1 TOft— - TOft - — — 



(3.6) 



It is important to note that even without viscous dissipation, the orbital decay naturally provides 
a dynamical tidal lag. This is because the star does not respond instantaneously to the changing tidal 
potential. For ui a » ft, this lag angle is 



Aft 

u 2 a t D ' 
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(3.7) 



(recall ASl ~ 0, since the spin is assumed to be small). This should be compared with the standard 
viscosity-induced tidal lag (cf. eq. [1.2]) 

Ail 

Ovisc ~ — 2T • (3-8) 

For most likely viscosities (see §8), we have to << t v i sc , thus the viscous tidal lag is negligible compared to 
the dynamical lag induced by the changing tidal potential. 

3.2 Resonance Amplitude 

As the orbital separation decreases, b a (t) increases according to equation (3.6). We can estimate how 
close to resonance this equation can still be valid, from which we can obtain a scaling relation for the 
maximum amplitude during a resonance. 

Suppose equation (3.6) starts to break down when 2fi ~ (1 — e)u> a . Note that e is not determined by 
uj 2 a — 4fi 2 ~ Q, which gives e ~ t or b/tD (evaluated at the resonance radius), where t or b — 2ir/tt is the orbital 
period; rather, equation (3.6) breaks down at a larger orbital radius, when the second term (the absolute 
value) in equation (2.17) becomes comparable to the third term, or, when the real part and the imaginary 
part in equation (3.6) becomes comparable, i.e., 

4m 3 Sl 2 Sl 2 22 

~ u) a - to Si . (3.9) 



ui 2 — m 2 il 2 

Therefore, the critical "closeness" to the resonance, before which equation (3.6) can be applied, is given by 



o\ 1/2 // x 1 / 2 



(3.10) 



where the subscript "a" implies values at the resonance radius D a . Evaluating equation (3.6) at 2Sl 
(1 — e)oj a , we obtain the amplitude at this near-resonance radius 



1 1 /0 1/2 



(3.11) 



The maximum amplitude |6 Qimaa; | during a resonance is of the same order. From equation (3.2) for the 
equilibrium tide, we have 

1/2 

(3.12) 







1 


(-) 




OL 


e 


\ torb J 



The maximum amplitude attained during a resonance can also be calculated as follows. Before the 
resonance is reached, a a (t) is essentially forced to oscillate at frequency 2Afi (cf. eqs. [2.10], [2.14], [3.6]). 
But after the resonance, we expect a a (t) to oscillate at its own eigenfrequency u a . Thus, for post-resonance 
oscillation, instead of equation (2.16), it is more appropriate to define a function c a (t) such such 

a a (t) - {GM'W lm Q nl )c a {t)e-^^\ (3.13) 



where s m = 1 for m = +2 and s m = — 1 for m = —2. Thus we have 

b a (t) = c a (t)e lm ^ lSMt . (3.14) 
Substitution of (3.13) into equation (2.10) yields 

c a - 2is m u a c a = J_ e --*+- m u^ (3 15) 

Integrating (3.15) with respect to time, and neglecting the c a term, we obtain the post-resonance tidal 
amplitude 

When tu >> l/w a , we can approximately set the upper (lower) limit of the integral to oo (— oo), and the 
integral can be evaluated using the stationary phase approximation. The maximum amplitude is then given 

by 

1 / \ 1/2 

' ; ' (m = ±2). (3.17) 



This has the same scaling as equation (3.11). Using equations (2.12), (2.15), (3.3) and (3.17), we have 



1 /57T 



16 



5/6 



1 /57r\ 1/2 /D \ 5/4 / Rr 2 \ 5/4 / 2 \ 5/4 

1 (bTT\ ( U a \ ( He \ q -i/2(_J_} (3.1 8 ) 



2 21 /4 ^ 6 J \RJ \GM J H \l + q 

5/4 



= 5.42 (D aA0 R 10 Mr.l) 5/i q~ 1/2 U 



where D Q) io = D a /(WR), R m = R/(10 km), M 1A = M/(1.4M ), and M Q is the solar mass. We see that 
the dimcnsionless resonance amplitude is larger for larger values of D a /R and Rc 2 /(GM). This is because 
the orbital decay timescale is larger for larger Dc 2 /(GM) (cf. eq.[2.13]), and therefore the binary spends a 
longer time at a particular resonance radius. For neutron stars, the ratio Rc 2 / (GM) ranges from 4 to 8 (e.g., 
Arnett & Bowers 1977); for white dwarfs, Rc 2 /(GM) ~ 10 3 - 10 4 . 

3.3 Numerical Solutions 

The dimensionless tidal amplitude b a (t) can be determined more accurately by numerical intergration. 
Define 

b a (t)=b^(t)+ib^(t), (3-19) 
where b& and b^ are real. Equation (2.17) then gives 

&£> + 2mO&« + (u, 2 m 2 tf)b^ + mf2&« = 

&« - 2mn&w - m n&w + ( w 2 - ™ 2 r! 2 )6« = o. 
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Numerical integration can be started at a distance larger than the resonance radius, where {D—D a )/D a » e 
(cf. eq.[3.10]). The initial condition for and bQ can be obtained from equation (3.6), i.e., 



D 



b (r) = 
"a 

&w = 

6« = t 

b^ ~ 0. 



2m 2 nn 



b (r) 



(3.21) 



Note that solutions for m — 2 and m = — 2 are related by 



b (r) _Ar) Ai) 
u m=2 ~ u m=-2> u m = 



m=2 



u m=-2- 



(3.22) 



Thus both m = — 2 and m — 2 modes are equally excited. 

In Figure 1 we show a typical example of the numerical result. We see that for larger orbital separation, 
D > D a , the tidal amplitude indeed evolves according a a (t) oc e ^ 4m * (cf. eqs.[2.16], [3.6]). Past the 
resonance, the tidal amplitude approximately evolves according to a a (t) oc e ^ ls ^ct _ rp ne absolute amplitude 
|c Q (t)| or |6 a (t)| does not exactly stay at a constant value after the mode is excited. This is because the 
tidal force continues to act upon the oscillation mode, although the net energy transfer is negligible after 
the resonance. Indeed, from our numerical results, we find that the average post-resonance amplitude |c a (t)| 
agree very well with equation (3.18), except when D a /R <> 3 — 4. 

Also note that the maximum amplitude given by equation (3.18) is not attained instantaneously at 
D = D a . From equation (3.16), the "duration of the resonance" 5t a , during which the resonance is effective, 
is given by 



St a 



jdte 



-im<$>-\-is m uj a t 



1/2 



1 

-^t or bto 



1/2 



(3.23) 



(see also RG2). The corresponding change in the orbital separation 5D a is 

1/2 



D a 



~ 0.05L» o 



l + q 



1/2 



{D a , w R^M lA ) 



5/2 



1/2 



(3.24) 



This expression agrees with the numerical result shown in Figure 1. 

4. EQUATION OF STATE AND CONVECTIVE STABILITY 
OF NEUTRON STARS 

Having understood the dynamical aspect of the problem of tidal excitations in §3, we now proceed to 
discuss the properties of the g-mode oscillations of neutron stars. 
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The restoring force for g-mode oscillations in the core of a cold neutron star is the buoyancy due to 
the chemical composition gradient, as first identified by RG1. However, the calculations of RG1 were based 
on an approximate (and inconsistent) ansatz for the Brunt- Vaisala frequency. In this section, we discuss 
how g-mode oscillations depend on the properties of nuclear matter, and how the Brunt- Vaisala frequency 
can be obtained self-consistcntly from some recent microscopic equations of state (EOS). Our calculations 
of g-modes will be presented in §5. 

4.1 Composition Gradient and Convective Stability 

The existence of stellar g-modes is closely related to the convective stability (stable stratification) of the 
star. For a sufficiently old neutron star, we can assume zero temperature and zero entropy throughout the 
star. Thus buoyancy can only arise from a composition gradient. Consider a blob of matter in equilibrium 
with its surrounding, with pressure P, the number of electrons per baryon Y e , and density p(P, Y e ). Imagine 
the blob is displaced upwards (against the gravity) adiabatically by a distance dr, where the surrounding 
matter has pressure P', composition F e ', and density p(P', V e '). The blob is in pressure equilibrium with the 
surroundings, but its composition is still Y e during the adiabatic process, since the timescale for the weak 
interaction, which is responsible for the change of Y e , is much longer than the dynamical timescale. The 
displaced blob density is therefore p(P',Y e ). Convective stability requires 

p(P',Y e )>p(P'X), (4.1) 

which gives the stability criterion (the Ledoux criterion) 

dp \ ( dY t 



8Y e J p \ dr 



< 0. (4.2) 
The gravity mode oscillation frequencies are closely related to the Brunt- Vaisala frequency N, defined 



as 

T2 „2'1 1 



=9 Z b-3 > (4-3) 



(e.g., Cox 1980). Here g = |g| is the local gravitational acceleration, c s is the adiabatic sound speed, as given 

(the subscript "s" means "adiabatic", which in this case implies constant composition), and c e is given by 

2 dP 

(the subscript "e" stands for "equilibrium"). We can easily rewrite the Brunt- Vaisala frequency (eq. [4.3]) 
as 

'(£),(£)--? &),(£)■ 

n 



where we have used the hydrostatic equilibrium condition dP/dr = —pg. Thus we see that the stability 
criterion (4.2) is equivalent to TV 2 > 0, or 

/ dP \ { dY \ 

c 2 - c 2 e = - —— > 0, (convective stability). (4.7) 



9Y e J p \dp 

We can now easily understand the existence of crustal g-modes (Finn 1987) and core g-modes (RG1) 
in a cold neutron star. In the outer crust, the matter consists of an electron gas embedded in a lattice of 
nuclei. As p increases, the nuclei become more neutron-rich (Baym, Pethick & Sutherland 1971), so that 
dY e /dp < 2 . On the other hand, since the pressure mainly comes from the electron gas, for a given total 
density, the pressure increases as the electron fraction increases, i.e., (8P/dY e ) p > 0. Therefore we have 
c 2 — c 2 > 0, and crustal g-modes are implied. 

In the core of a neutron star, the signs of dY e /dp and (dP/dY e ) p are opposite to those in the outer crust. 
Consider a simple free npe (neutron, proton and electron) gas model for the core (e.g., Shapiro & Teukolsky 
1983, RG1). The electron fraction Y e increases with density, dY e /dp > 0. But since the pressure is mainly 
provided by the neutrons, we have (dP/dY e ) p < 0. Again we see N 2 > 0. More detailed discussions of these 
points based on general EOS of nuclear matter are given in §4.2. 

4.2 Nuclear Equation of State and Brunt- Vaisala Frequency 

We see from §4.1 that calculations of the core g-modes of a neutron star requires knowledge of both c s 
and c e , in addition to the density-pressure relation. Most of the early EOS's (e.g., Arnett & Bowers 1977) do 
not provide enough information to calculate the sound speed and to determine the Brunt- Vaisala frequency. 
Here we discuss how to obtain these information from a general parametrization of nuclear matter. Such 
paramctrization has been used in several recent microscopic nuclear EOS's. 

Consider the high density liquid core of a neutron star consisting of neutrons, protons and electrons in 
(3— equilibrium. In general, as a function of the baryon number density n and the proton fraction x = n p /n, 
the energy per baryon of the nuclear matter can be written as (Lagaris & Pandharipande 1981; Prakash, 
Ainsworth & Lattimer 1988; Wiringa, Fiks & Fabrocini 1988, hereafter WFF) 

E n (n,x) =T n (n,x) + Vo(n) + V 2 (n)(l-2x) 2 , (4.8) 

where 

T n (n,x) = |^_(3^ 2 „)V3 [a; 5/3 + (1 _ x)5 / % (49) 
5 2m n 

is the Fermi kinetic energy of the nucleons, m n is the nucleon mass. The last two terms in equation (4.8) 
represent contributions due to nucleon interactions: Vq mainly specifies the bulk compressibility of the matter 
(hence the pressure-density relation) , and V2 (the "symmetry potential" ) is related to the symmetry energy 



2 Note that the change of Y e in the crust is actually discontinuous. 
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of nuclear matter 3 , which plays an important role in determining the equilibrium proton fraction. The 
parametrization of dense nuclear matter by equation (4.8) is expected to be generic, and it can simulate 
many recent microscopic calculations. 

The energy per baryon of relativistic electrons is given by 



T e {n,x e ) = lhc(37r 2 n) 1/3 xi /3 , 



(4.10) 



where x e = n e /n = Y e . Charge neutrality requires x — x e . The total energy per baryon in the liquid core is 
then 

E(n,x) = E n (n,x)+T e (n,x e ). (4.11) 
Minimization of E(n,x) with respect to x, i.e., dE/dx = 0, yields the condition for /3-equilibrium 

M« - M P = Me, (4.12) 

where jtij's are the chemical potentials of the three species of particles, with 

QE H 2 

and 

p e = &c(37r 2 n) 1/3 *e /3 - 



^ = 4(1 - 2x)V 2 + -^(37r 2 n) 2 / 3 (1 - x) 2 ' 3 - r 



2/3 



(4.13) 



(4.14) 



The equilibrium proton fraction x(n) — x e (n) can then be obtained by solving equations (4.12)-(4.14). The 
mass density and pressure are determined via 

p{n,x) = n[m n + E(n,x)/c 2 }, 

\ 9 dE(n,x) 2r* 7i , r . . . „ (4-15) 

P(n, x) = n 2 ^ > = —T n + -T e + n 2 [V£ + F 2 '(l - 2a;) 2 ] , 

where Vq = dVo/dn, etc. The equilibrium EOS is then obtained when p and P are evaluated at the 
equilibrium composition x = x(n). The adiabatic sound speed c s is given by 
9 dP n dP 

(4.16) 



dp p + P/c 2 dn 

{ Y Tn + \ Te + 2n + V2 (1 ~ 2x)2] + n ' + V2 " {1 2x)2] } ' 



p + P/c 2 

The difference between c 2 and c 2 is given by 



2 2 
C„ - C = 



(dP_ dP\ _ 



/dP\ dx 



p + P/c 2 \dn dn J p + P/c 2 \dx J dn 



p + P/c 2 



d , 

— {fj, e + p p - p n ) 



dx 
dn 



(4.17) 



3 The symmetry energy E s is defined such that E n (n, x) = E n (n, 1/2) + E s (n)(l — 2x) 2 H — •. Therefore, 
E s {n) = 5T n (n, l/2)/9 + V^ra). The experimental value is E s ~ 30 MeV at the saturation density n s = 0.16 
fm- 3 . 
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From equation (4.12), we have 



dx 
dn 

Equation (4.17) then becomes 



2 2 

ct-ci = 



p + P/c 



1 2 



9n 



(He + [i p - //„) 



9x 



Clearly, convective stability (cf. eq.[4.7]) requires d(fj, e + [i v — fj, n )/dx > 0, i.e., 



^Me + 8V 2 + ^( 3 7r 2 n) 2 / 3 [(1 - x)-V* + x- 
6x 6 2m n L 



1/3 



> 0. 



(4.18) 



(4.19) 



(4.20) 



Thus unless V 2 is extremely negative, core g-modes always exist. On the other hand, zero-temperature 
neutron stars must always be convectively stable, since for the cold nuclear matter in the ground state, there 
is no energy available to sustain convective motion 4 . Therefore equation (4.20) provides a constraint on 
the nuclear equilibrium state. 

Finally, we note that just above the nuclear density, when fi e exceeds the muon rest mass m M , muons 
are energetically allowed via n <-> p + p,. The /3-cquilibrium condition becomes 



Hn - Mp = Me = M/i = [ TO u c4 + ft 2 C 2 (37T 2 nX M ) 2/3 ] 1/2 , 



(4.21) 



with The expressions derived in this section can similarly be modified. However, there are 

greater uncertainties in our knowledge of nuclear EOS, especially T^(n). For definitcness, we will not include 
muons in our calculations. 

4.3 Neutron Star Core EOS Models 

In this paper, four nuclear EOS's will be considered. The first three are based on recent microscopic 
calculations of WFF. For n > 0.07 fm~ 3 , the functions Vb(n) and V^(n) obtained from different nucleon 
Hamiltonians have been tabulated in Table IV of WFF. However, to facilitate smooth derivatives of Vq and 
V2 with respect to n (which are needed to determine c s and especially c 2 — c 2 ), we find it more convenient 
to use approximate fitting formulae, which we give below: 

(1) Model AU: This is the EOS based on nuclear potential AV14+UVII of WFF. V and V 2 (in MeV) 
are fitted as 

V = -43 + 330 (n-0.34) 2 , 

(4.22) 

y 2 = 21n ' 25 , 
where n is the baryon number density in fm~ 3 . 

(2) Model UU: This is based on potential UV14+UVII of WFF. The fitting formulae arc (in the same 
units as in eq.[4.22]) 

V = -40 + 400(n- 0.3) 2 , 

(4 - 23) 

V 2 =42n ' 55 . 



4 I thank Ed Salpeter and Andreas Reisenegger for pointing this out to me. 
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(3) Model UT: This is based on potential UV14+TNI of WFF 5 . The fitting formulae are (in the same 
units as in eq.[4.22]) 

V Q = -42 + 350 (n-0.28) 2 , 

(4.24) 

V 2 = 18-130 (n-0.29) 2 . 

Once the functions Vo(n) and ^(n) are given, all the relevant EOS parameters needed for the calculations 
of g-modes can be obtained using the expressions in §4.2. Some of the results are shown in Figure 2. 

The above fitting formulae are valid only for n <y 1 fm~ 3 . This is adequate for our purpose, since we 
will be mainly concerned with M = 1.4M Q NS's, for which all three EOS's yield a central density that is 
less than 1 fm -3 . For higher density, both AU and UU violate causality, indicating considerable theoretical 
uncertainties. We note that for AU, the fitted function V 2 does not reflect the exact behavior of V2 in Table 
IV of WFF (e.g., V2 in that table shows non-monotonic behavior for 0.07 fm~ 3 < n < 0.7 fm~ 3 , while our 
fitted formula is monotonic in n). However, there is much greater uncertainty in the three-body nucleon 
interaction potential (especially its isospin dependence) that determines V2{n) itself. Also note that for 
model UT, the quantity c 2 — c 2 becomes negative at very high density. This implies that the model is not 
reliable in such a high density regime (see discussion following eq. [4.20]). 

(4) Model UU2: Our fourth EOS has the same P(n), p(n), and c e (n) as UU. However, we shall adopt 
the ansatz similar to that of RG1 for c 2 — c 2 in order to determine to what extent the g-mode properties are 
affected. This ansatz is based on the following consideration. 

For free npe gas (i.e., set Vo = V 2 = for the equations in §4.2), /^-equilibrium gives (see eqs.[4.12]-[4.14]) 

x = 3 ^ n ( _A \ [(1 _ )2 /3 _ x 2/3 ]3 „ ^2 ( 



5.6 x 10" 3 - 



(4.25) 



2.8 x 10 14 g/cm 3 ' 

where the second equality follows from x « 1. The sound speed is given by (eq.[4.16] with Vq = V 2 = 0) 



2 1 dP 2 H 2 



m n dn 3m. 
From equations (4.13)-(4.14) and (4.17), we have 



— 2^< 3 * 2 ») 2/3 ( 1 -f)' <« 6 » 



^2 c 2 1 dx 

m n dn 



n 2n 



(4.27) 



Using equations (4.12), (4.26)-(4.27), we then have 



This agrees with the expressions given in RG1 (see their eqs. [28]-[29]). 



5 The results tabulated in Table IV of WFF contain slight errors (Wiringa 1993, private communication). 
We have used the updated results provided by Bob Wiringa. 
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Now the ansatz for UU2 is to use equation (4.28) to obtain c s from c e and to use x(n) from equation 
(4.25), although both equations are only valid for free npe gas EOS model 6 . This ansatz is clearly not a 
consistent procedure to obtain c 2 — c 2 and hence the Brunt- Vaisala frequency. In a realistic description of 
nuclear matter, the sound speed c s is determined by the compressibility of the matter, while the quantity 
c 2 — Cg is mainly related to the symmetry energy. Therefore, equation (4.28) cannot hold in general. Indeed, 
from Figure 2 we see that the ratio (c 2 — c 2 e )/c 2 s in UU2 is qualitatively different from that in UU, although 
the other quantities are the same. Nevertheless, in view of the great uncertainties in our knowledge of T^(n), 
it is worth exploring the g-mode properties based on this ansatz. 



We now discuss our calculations of the g-mode oscillations of a isolated cold neutron star based on the 
EOS models of §4. The resulting mode frequencies ui a and the tidal coupling coefficients Q n i (eq. [2.11]) will 
be used in §6 to determine the energy transfer and angular momentum transfer due to resonant excitations 
in coalescing neutron star binaries. 



In this paper, we use Newtonian linearized equations to calculate the oscillation modes. The use of 
Newtonian equations is consistent with our Newtonian description of tidal interactions (§2). For f-mode, 
general relativistic effects are expected to modify our results of oscillation frequencies by not more than 
GM I (Re 2 ) ~ 20%. Linear oscillations of nonrotating stars in Newtonian theory have been discussed exten- 
sively in the literature (e.g., Ledoux 1974, Cox 1980). Here, for completeness and to clarify the notations, 
we give some of the relevant equations below. 

The basic equations that govern the oscillations of stars are the equations of hydrodynamics, i.e., the 
continuity equation, the equation of motion, energy equation (which specifies the adiabatic nature of the 
oscillations), and the Poisson equation for the gravitational potential. The linearized equations can be written 
in terms of the radial component of the fluid displacement £ r (cf. eq.[2.8]), the Eulerian perturbations 
of density, pressure and gravitational potential, Sp, 5P and S(j). Separating out the angular dependence 
Yi m (9, <f>), these linearized equations are (here and throughout this section, the mode index a = {n, I, m} is 
suppressed) 



6 RG1 have also considered the ansatz using more realistic x(n) in equation (4.28). This is also not a 
consistent procedure. 



5. CORE g-MODE OSCILLATIONS AND TIDAL 



COUPLING COEFFICIENTS 



5.1 Basic Equations 




(5.1) 
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where u is the eigenfrequency of the mode, and 



r, = ^, A-—, t^^y-y-^n. (5.2) 

Here, c 2 , N 2 can be calculated using equations (4.3), (4.16)-(4.17). The transverse component of the dis- 
placement vector (see eq. [2.8]) is given by 



— — + 5(j){r) 



(5.3) 



The eigenmodes can be determined by solving equation (5.1) with appropriate boundary conditions. 
For r — > 0, regularity requires 

r 2 f = -^(Y + <S> )r l+1 , — = Y r l , 5<t> = 3> r l , (r -> 0), (5.4) 

uj z p 

where Y Q and $ are constants. At the stellar surface, the Lagrangian perturbation of the pressure, AP, 
must vanish, i.e., 

AP = 5P — pg£ r = 0, (r = R). (5.5) 
Also, the gravitational potential and force must be continuous across the surface. This gives 

O 7 I 1 

—54* + ——(50 = -47rGp£ r , (r = R). (5.6) 
or r 

Numerically, we fix Y = 1, which simply sets a scale for the eigenfunctions. Equations (5.1) are 
integrated outward from r — > with boundary conditions (5.4), using a fifth-order Rungc-Kutta scheme. 
The outer boundary conditions at r = R (eqs.[5.5]-[5.6]) can be matched using a standard shooting routine 
(Press ct al. 1992), therefore the eigenvalues lo 2 and $ Q can be obtained after an initial guess. The number 
of nodes in the functions is counted to determine the orders of the modes. After to 2 is obtained for a specific 
mode, the eigenfunctions are normalized according to equation (2.9), i.e., 

i-R 

pr 2 dr [C{rf + 1(1 + l)^(r) 2 ] = 1. (5.7) 

Then the tidal coupling coefficient Q n i is calculated using equation (2.11). To test our numerical code, we 
have done calculations for polytropic stellar models, with P = Kp r , where the polytropic index T and the 
adiabatic index Ti > T are constant throughout the star. For a few specific values of T and Ti, the mode 
frequencies and Q n i have already been tabulated in the literature (e.g., Lee & Ostriker 1986). Our results 
for the polytropic models agree with these tabulated values. 

We construct equilibrium neutron star models by solving Newtonian hydrostatic equations. We use 
Newtonian equations for consistency with our use of Newtonian equations in calculating oscillation modes 
and in describing the tidal interaction. The EOS's we have used are those discussed in §4.2-4.3 for n > 0.07 
fm~ 3 , where the pressure, sound speed, and c 2 — c 2 (hence the Brunt- Vaisala frequency) are all calculated 
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self-consistently (except model UU2). For lower density crustal regions, we use the Baym, Bethe & Pethick 
(1971) EOS above neutron drip (p^4x 10 11 g/cm 3 ), and then join onto the Baym, Pethick & Sutherland 
(1971) EOS below neutron drip. In all of our calculations, we set c 2 s = (? e in the crustal region, therefore 
effectively suppressing the crustal g-modes while concentrating on the core g-modes. Since the density in 
the crustal region is much smaller, from equation (2.11) we see that the tidal coupling coefficients for these 
crustal modes are likely to be very small (as found by Shibata 1993, who used highly simplified neutron star 
models). 



Table 1 gives our numerical results for the oscillation frequencies and tidal coupling coefficients Q n i of 
the f-mode and the first two g-modes for neutron stars based on the four EOS models discussed in §4.3. 
Note that since the gravitational wave frequency is twice that of the orbital frequency f2, and resonances 
occur when uj a = 20, the listed mode frequencies also correspond to the gravitational wave frequencies at 
resonances. We see that the four neutron star models considered have very similar stellar radii for the given 
canonical mass M = 1.4M Q , and their f-mode properties are also very similar. This is because the four 
EOS's of §4.3 predict very similar bulk properties (e.g., pressure, compressibility) for the nuclear matter. 
However, the g-mode properties are very different for the four models. In particular, although UU and UU2 
yield the identical neutron star mass-radius relations, the g-mode frequencies predicted by them differ by 
large factors. These differences reflect the sensitive dependence of g-modes on the symmetry energy of the 
nuclear matter. 

From Table 1 we see that the frequencies of g-modes are much smaller than that of f-mode. This is a 
direct consequence of the small difference between c s and c e (cf. Fig. 2). Also, \Q n i\ becomes smaller as the 
radial order of modes increases. This is because of the strong cancellation in the integral for Q n i (eq. [2.11]), 
as already noted by Cowling (1941). Indeed, using the equation of continuity, 



Since the order of a mode corresponds to the number of nodes in 8p n i {r) , we clearly see that for the higher- 
order mode \Q n i\ is smaller. Based on our numerical values of Q n i for the first few g-modes, we find that for 
all NS models considered, \Q n i\ decreases as uj n i decreases faster than \Q n i \ oc ui^f- 

The frequencies of the first ten quadrupole g-modes are given in Figure 3. They generally lie in the range 
of 10 Hz to 100 Hz. An approximate asymptotic expression for the frequencies can be easily obtained. For 
high-order g-modes (n >> 1), the perturbation of the gravitational potential 6<j) can be neglected ("Cowling 



5.2 Numerical Results and Discussions 



(5.8) 



the expression for Q n i (eq. [2.11]) can be rewritten as 




(5.9) 
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approximation"). Eliminating SP/p from the first two equations in (5.1), and noticing that lo — > 0, we have 

"ra^lf-)-^^). <«»> 



Or 2 ' uj 2 dr \ p ) uj 2 r 2 
Assuming r 2 £ r ~ e lkr , we obtain the WKB local dispersion relation for high-order g-modes 

fc W 2 - ^W) 2 - (5-11) 
The eigenvalues for w 2 are then determined by the stationarity condition 

k(r)dr = (n + C)ir, (5.12) 

o 

where C is a constant of order unity. Thus, asymptotically, the high-order g-mode angular frequencies are 
given by 

The values of Q n i for the higher-order g-modes are harder to calculate numerically as they become much 
smaller. Indeed, in evaluating the values of Q n i, it is necessary to know the eigenf unctions of the g-modes 
very precisely, since a small f-mode contamination can significantly change the obtained values (Reisenegger 
1994). To obtain meaningful numbers of Q n i for the g-modes, it is also important to implement self-consistent 
hydrodynamical equations in the numerical calculations, so that the Hermiticity of the oscillation operator 
C is maintained. For example, we have found that significant systematic errors in Q nl can be introduced 
if c 2 — dP/dp is not satisfied to a high accuracy when interpolating tabulated EOS. For the g-modes given 
in Table 1, we have checked that the orthogonality relation (2.9) is satisfied to sufficient accuracy, therefore 
our values of Q n i for these g-modes are not contaminated by the f-modc. 

We have also done calculations of the oscillation modes using the Cowling approximation. We find, as 
expected, that the g-mode frequencies are nearly unaffected. However, for a given stellar model, the Cowling 
approximation leads to a larger f-mode frequency (by ~ 50%) and larger values of \Q n i\ for the first few 
g-modes (by as large as 



6. TIDAL ENERGY AND ANGULAR MOMENTUM 

In this section, we combine our results of §3 and §5 to calculate the amplitude of distortion, energy 
transfer and angular momentum transfer to the neutron star due to tidal interaction. The effects on the 
orbital phase will be considered in §7. 

6.1 Tidal Distortion Amplitude 

A measure of the tidal distortion of the neutron star is the quantity a a (cf. eq.[2.6]). For equilibrium, 
quasi-static tide (m = 0), equation (3.2) gives 

iT 3 



|a a ,eq| ~ qiJ nl \W 2a Qnl\ [p] ■ C 6 - 1 ) 
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where to n i is the mode angular frequency in units of (GM/i? 3 ) 1 / 2 . The dominant tidal distortion comes from 
the f-mode, for which u n i ~ 1 and Q n i ~ 1 (see Table 1), thus 



h /i? x3 



Ke 9 | ~ =? (^J , (6.2) 

where /to is the typical height of the tidal bulge (see eq. [2.19]). 

Now consider the dynamical tides (m = ±2). Before the resonance, the distortion amplitude associated 
with a particular mode increases as the orbit decays (see eq. [3.6]). After the resonance, the amplitude stays 
nearly constant. From equations (2.16) and (3.18), the tidal distortion amplitude resulting from the resonant 
excitation is given by 

i i n -s/ein \( Rc2 \ 5,i W 2 X5/6 



-3 / fnl \ ( \Qnl\ \ M ~5/6 1/2 7 



(6-3) 

where /„; = ujni/2-K is the mode frequency which is also equal to the frequency of the gravitational wave 
when the resonance occurs. The reference values /„; = 100 Hz and \Q n i\ = 0.0003 have been chosen to 
approximately agree with the lowest order g-mode in model UU (see Table 1). Note that equation (6.3) is 
valid for both m = 2 and m = —2, and it is half of the result given by RG2 (see their eq. [14]) 7 . Since 
\a a \ « 1, the linear approximation we have adopted throughout the paper is valid. 

6.2 Tidal Energy Transfer 

Next we look at tidal energy transferred to the neutron star. The kinetic energy of the star is given by 

EW-lJt^.g-lZM, (6.4) 

a 

where we have used equations (2.6) and (2.9). Similarly, using equation (2.7), the potential energy (including 
the internal energy) of the oscillations is given by 

E P {t) = \ J<Px?-£-£=±'£ i u>l\a a \ 2 , (6.5) 



7 In the PT formalism used by RG2, the contributions from m = 2 and m = — 2 terms are not equal, while 
it is clear from our analysis (see eqs. [2.10], [3.16] or [3.22]) that the m — 2 and m = — 2 modes contribute 
to tidal distortion and energy equally. The difference comes from the fact that contributions from m = 2 
and m = — 2 modes have been re-grouped to derive the final expressions of PT. Thus, although explicit 
evaluations in RG2 are made only for I = m = 2, RG2's results actually include contributions from both the 
m — 2 and m = — 2 modes, since the other term in PT's formula are much smaller. 
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(e.g., Shapiro & Teukolsky 1983, pl41-143). The total tidal energy of the neutron star is E = E k + E p . 
Substituting equations (2.16) and (3.13), and using the natural units introduced at the end of §2, we have 



; a | > 



= \ ^{M'WimQnifK - imnb a \ 2 = 1 ^{M'W lm Q nl f\c a - is m uj a c 

a a (6.6) 

a a 

Another way to consider the tidal energy of the neutron star is to look at the rate of energy transfer: 

E = - Jd 3 xp^-VU*. (6.7) 
Using equations (2.2), (2.6), (2.11) and (2.16), we have 

E = J2(M'W lm Q nl )^e im9 a a = ^{M'W^Q^) 2 ^{b a - imSlb a ). (6.8) 

a " a 

It is easy to check that E in this expression is real (as it should be) , since the contributions to the imaginary 
part from the m = +2 term and that from the m = — 2 term cancel. Moreover, using equation (2.10), we 
can show that this expression for E is consistent with equations (6.4)-(6.5) for Ek and E p . 

For equilibrium, quasi-static tide (m = 0), when orbital decay time to » l/u>nz, the kinetic tidal energy 
is small compared to the potential energy (thus "quasi-static" ) . The equilibrium tidal energy associated with 
a particular mode is then 

1 GM 2 2 2/ x2 fR\ 6 , , 

E ^=2^ q ^ {W2oQnl) \D) ■ (6 ' 9) 

Again, the dominant equilibrium tidal energy comes from the f-mode 

GM 2 , /R\ 6 

E ^~—<i\n) . ( 6 - 10 ) 

which is the expected result (e.g., Lai et al 1993). 

Now consider the energy transfer AE n i due to the resonant excitation of a specific g-mode. Since c a (t) is 
approximately constant after the resonance, the kinetic energy and the potential energy are equal. Equations 
(3.18) and (6.6) then yield 

n 2 ( GM 2 \ i/3 n 2 fRc 2 \ 5/2 J 2 \ 5/3 ( 6 - n ) 



512 {-fT)^ 1 Qnl \GM) q \~TTq 



Note that the factor of 2 in the first equality comes about because the m = 2 and m — — 2 terms contribute 
equally to the excitation energy, and we have included both of them in (6.11). The ratio of AE n i to the 
orbital energy 

GMM' ,„ 
Eorb^ 25", (6.12) 
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is given by 

\E or b\ 128 



AE nl tt 2 _ 1/3 ^ 2 /i?c 2 \ 5/2 ' ° X4/3 



(6.13) 

1/3 / , n I x 2 / x 4/3 v ; 



,1 + 9, 

where |-E rb| is evaluated at D = D n \. This ratio plays an important role in determining the effect of 
resonant energy transfer on the orbital decay rate (see §7). Since AE n i « \E or b\, the resonance has only 
a small effect on the orbit, and the use of the trajectory of ideal two point masses (cf. eq.[2.12]) is a good 
approximation. Our result for AE n i agree with that of RG2 (see their eqs. [15] and [17]), but note that our 
equations (6.11) and (6.13) include both m = 2 and m = — 2 modes (see footnote 7). 

6.3 Angular Momentum Transfer and Dynamical Spin-up of the Neutron Star 

Associated with the tidal energy transfer is an angular momentum transfer to the neutron star. This 
is possible because of the dynamical tidal lag induced by the decaying orbit, as already noted in §3.1 (see 
eqs. [3.6]-[3.7]). The z— component of the tidal torque on the neutron star per unit mass is given by 

T z = -e z -(rxV)U = - — U. (6.14) 

Thus the total torque is 

N z = Jd 3 x(p + 5p)r* = Jd 3 x5pr*. (6.15) 
Using the continuity equation (5.8), we have 



= 2(M'Wi m Q nI ) 2 5 J jT (-im& a ) - J2( M ' W i™Q^ 2 W+T mb « 



(6.16) 



Note that this is non-zero since the contribution from the m = +2 term and that from the m = —2 term add 
up (cf. eq.[3.22]). Clearly, the equilibrium tide (m = 0) does not contribute to angular momentum transfer. 

Prior to the resonance, the tidal torque can be evaluated using equation (3.6). In order of magnitude, 
the torque associated with a specific mode is 

K, a ~^f q 2 ^Ql t (6.17) 

where the dynamical tidal lag angle S a is given by equation (3.7). The dominant f-mode contribution is 

GM 2 2 /i?\ 6 r GM 2 2 /R\ e n , C10 , 

- — q [d ) — q [d) ^r D (6 - 18) 

It is easy to see that this non-resonant dynamical tidal torque does not appreciably affect the neutron star 
spin. In particular, this torque is not sufficient to synchronize the spin with the orbital angular velocity: the 
timescale for synchronization is given by 

tsyn ~ ~A4~o \R ) tD > (6 - 19) 

22 



which is much larger than the orbital decay time to unless D becomes comparable to R. 

Now consider the total angular momentum transferred to the star, AJ a , due to a particular resonance. 
From equation (6.16), this is given by 

AJ a = (M'W lm Q nl ) 2 J dt^(-imb a ), (6.20) 

where only the real part should be included (recall that the imaginary part from the m = 2 term and that 
from the m = — 2 term cancel). Using equations (3.14)-(3.15), we have 



AJ a = {M'W lm Q nl ) 2 {-im) J dt-±_ 

= (M'WimQni) 2 (-im) J dtc a (c* a + 2iu a s m c* a ). 



im& — iu> a s m t 

(6.21) 



Integrate by parts, and notice that the term proportional to c a c* can be neglected, since after the resonance, 
c a is nearly constant. Adding up contributions from both m — 2 and m — — 2 terms, the total angular 
momentum transfer A J n i due to the resonance is given by 

AJ n i ~ 2(M'Wi m Q n i) 2 uj n i\m\\c a . max \ 2 = AuJ n i\a a , ma x\ 2 ■ (6.22) 

Using equation (6.3), we have 

2 / R 2 \ 5 / 2 / n \ 5/3 



Note that the energy transfer and angular momentum transfer are related by 

AJnl = ^L = ^L, (6 . 24) 

where VL n i = to n i/2 is the orbital angular velocity at the resonance. 

The effect of this resonant angular momentum transfer on the spin of the neutron star can be considered 
if we assume that the angular momentum goes into uniform rotation of the star. The resulting change of 
the spin rate due to the resonance is then given by 

Att s ,„ ; = (6.25) 

where I — 2kMR 2 /5 is the moment of inertia of the star, k < 1 is of order unity. The ratio of AQ s _ n i and 

n n i is 

-5/3.2 (RC 2 \ 5/2 ( 2 \ 5 / 3 



Q n i 256k 



1.5 x lO" 4 /*- 1 [ I [ I M 



fnl \ ( \Qnl\ V ,,-5/3 „ ( 2 



(6.26) 



,100 UzJ V°- 0003 / V 1 + 9, 

Therefore, spin-up due to resonant angular momentum transfer is negligible. In reality, the angular momen- 
tum transferred does not manifest itself as uniform spin of the star, since the viscous timescale to achieve 

23 



uniform rotation is long. Thus the actual change of spin rate is even smaller than that given by the above 
expression. 



7. EFFECTS ON ORBITAL DECAY RATE: PHASE ERROR 

We now calculate the orbital phase error due to the resonant energy transfer and angular momentum 
transfer from the binary orbit to the neutron star. This is important because, as mentioned in the intro- 
duction, a small phase error in the theoretical gravitational wave template can destroy a possible detection 
using the matched filter technique. 

Because of the relation (6.24), a circular orbit will remain circular after the resonance. Therefore it is 
sufficient to consider the energy of the system only. Excluding a constant, the total energy of the system, 
E tot , can be written as a sum of the orbital energy E orb and the stellar energy E, i.e., E tot = E orb + E. Thus 

• _ dE tot ■ f dE orb dE\ ■ 

Etot ~ ~dD D \~W + dD) D - (1A) 
The number of orbital cycles, N orb , can be obtained from 

Q dD n ,^fdE orb dE\ ,. 

The first term in (7.2) is simply the point mass result. The second term represents the error induced by the 
change of stellar energy: 

d(AN orb ) = ^-dD (^j /Etot- (7-3) 

This is a general expression, which applies to dynamical tidal effects as well as equilibrium tidal effects. Now 
consider the total accumulated error on the orbital cycles AN orbiTl i due to resonant excitation of a particular 
mode with angular frequency w„;. Since dE/dD is large only near the resonance radius D n i (see eq.[3.3]), 
we can integrate equation (7.3) to obtain 

a A , Qnl AE n i to AE n l 

AN orb , nl ~—-^- = - — —- (7.4) 

Z7r Etot '-orb \&orb\ 

where to = \E orb / Etot\ = \D/D\ (see eq. [2.13]), and t orb = 2tv/Q, is the orbital period. Note that all the 
quantities in equation (7.4) are evaluated at the resonance radius D n i. Here the negative sign implies that 
the effect of energy transfer is to make the binary coalesce slightly faster (i.e., slightly larger \D\) compared 
with the ideal point mass case. Substituting equations (2.13) and (6.13) into (7.4), we have 

„2N 5 



AN orb , n i * ~^-jQlt (|^) « (-^) 



(7.5) 

~ -4.3 x 1(T 4 I -rpT^j- ) 1^1 M-^ / 



-4 1 fnl \ ( \Q 



^lOOHzy VO-0003/ 14 w ^\l + q 
Note that this phase error does not accumulate instantaneously at D = D a (see §3.3). From equation (3.23), 
the number of orbital cycles that unfold during the period of effective resonance is given by 

~ ^ * 12 - {mkT" M ""^' 2 (ih)""- (7 - 6) 
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The small size of AN or b, n i in equation (7.5) implies that the phase errors due to resonant excitations of g- 
modes are negligible for constructing gravitational wave templates. Therefore, to a very good approximation, 
binary neutron stars can be treated as point masses during the inspiraling phase. Hydrodynamical tidal 
effects are dominated by the f-mode distortion, and they are important only when the neutron stars are close 
to contact (see eqs. [6.2], [6.10]). 

8. TIDAL HEATING OF BINARY NEUTRON STARS 

Having determined in §6 the tidal energy in the decaying binary neutron stars, a natural question 
concerns where this energy goes. Without damping mechanisms, this energy is stored in the form of potential 
energy and kinetic energy of various oscillation modes. There are basically three processes that can lead 
to damping of g-modes (RG1): gravitational radiation, relaxation toward chemical equilibrium by neutrino 
emission, and viscous dissipation. The damping timescales for all of these mechanisms are likely to be much 
larger than the orbital decay time of interest here, i.e., the last few minutes prior to the binary merger. 
However, viscous dissipation leads to heating of the neutron stars, thus changing the physical properties of 
the stars even before merging takes place. Therefore it is of interest of consider the heating processes in 
some detail. 

The standard estimation of viscous tidal dissipation goes as follows. The rate of energy dissipation is 

(8.1) 

where v is the kinematic shear viscosity (in units of cm 2 /s), and dvi/dxk is the strain tensor. The typical 
tidal bulge height is ho, and the bulge rises and falls at angular velocity 2 Ail (see eq. [1.2]). If the tidal 
distortion of the star varies over a length-scale comparable to the stellar radius, then the shear strain tensor 
is of order dvi/dxk ~ Ailho/R. Thus the viscous dissipation rate is 

E msc ~ vM (aA^) 2 = M{h ^ ? = vM{AQ) 2 q 2 (J) 6 . (8.2) 

This is simply the tidal kinetic energy Ek ~ M(h An) 2 divided by the viscous time t visc ~ R 2 /v. For a 
synchronized system (Afi = 0), there is no viscous dissipation. 

It is clear from the above derivation that this simple relation is mainly for the f-mode tide, i.e., higher- 
order modes and resonances have been neglected. However, higher-order mode oscillations involve more shear 
motion of the fluid, thus the length-scale and timescale for viscous dissipation are smaller for these modes. 
More importantly, when g-mode resonances occur, the viscous dissipation can start at a larger separation, 
at which the orbit decay is still slow. As a result, more heat can be generated during the whole inspiral. By 
contrast, dissipation due to the non-resonant f-mode oscillation is mainly effective for small orbital separation 
(see eq. [8.2]), at which binary decay is fast. 

Of course, the amount of viscous dissipation ultimately depends on the viscosity of the fluid, which is 
not at all well understood. Our main intention in this section is to consider the effect of g-mode resonances 
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E vis , ~ vM 



dvi 



dxi 



on the viscous dissipation rate. Although the stellar spin can be incorporated in our formalism without much 
difficulty, we shall focus on the non-spinning case, so that AQ = ft. In this case, equation (8.2) indicates 
that for non- resonant viscous dissipation E visc oc 1/D 9 . 

8.1 Viscous Damping Rate of Normal Oscillation Modes 

Here we consider the viscous damping rates of individual modes. To leading order, viscous stress does 
not change the oscillation frequency, but it induces mode damping. This damping effect is equivalent to 
adding an imaginary part to the oscillation frequency, i.e., uj' a = ui a + ij a . The mode energy damping rate 
is given by 

2 7q = ^pEL. (8.3) 

Here E a is the energy of the eigenmode: 

E a = 2E aM = <*%f d 3 xp C ■ & = ul ■ (8.4) 

where we have used the normalization (2.9). The viscous dissipation rate is given by (suppressing the mode 
index a) 

E visc = - d 3 xa ik v* k = - I d 3 x 



1 2 

2?7 K,fc + v ka - g<5ifcV • v| 2 + C| V • v| 



.5) 



(Landau & Lifshitz 1987) where Oi k is the viscous stress tensor 

Oik = f] {vifr + v k ,i - ^<5 ifc V ■ v^j + (5 ik \7 ■ v, (8.6) 

7] = pv is the dynamical shear viscosity, and ( is the bulk viscosity. 

Substituting v = — iu;^ and equation (2.8) into equation (8.5), the damping rate can be evaluated. In 
the absence of bulk viscosity, the result has been obtained by Higgins & Kopal (1968, see also Kopal 1978). 
For m = modes, 

lt h Jo r) = \[r*drr,F{r), (8.7) 



where 



*Xr)=2 +[^--1(1 + 1 



r 



dr 



(8.8) 



and 



For modes with m^O 8 , 



G(r) = -£- + -? -1(1 + 1)$-. (8-9) 
Or r r 



(shear) _ (HJH)I (shear) , R 

7 "*o - (z-H)! 7m =° ' (8 - 10) 



8 Our definition of Yi m is that of Jackson (1975), which is different from that used in Higgins & Kopal 
(1968). Therefore the results are different by a numerical factor. We have checked that for the f-mode of 
incompressible sphere (the Kelvin mode), our result using equations (8.7)-(8.9) agrees with the standard 
one, as first derived by H. Lamb; see Cutler & Lindblom (1987). 
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When bulk viscosity is present, we can add the bulk viscous dissipation to the damping rate, i.e., 

7 = -y( S ' lear ) + ryibulk) ^ (8.11) 

where 

7» = \I**<KGM\ 7 <™> = f^gf 7»- (8-12) 

To specify the dependence of 7 on the mode structure, independent of the actual values of viscosities, 
let us consider viscosities of the form 



v = m[ — ), C = 6v — • (8-13) 

\PnJ \PnJ 

where pat is some fiducial density. Such a density dependence is appropriate for the microscopic viscosity in 
the core of neutron star (§8.3). The damping coefficient can be written as (for the m — mode) 

(shear) _ VnJ^ ( M \ 2 y(shear) (bulk) _ C^R / M \ {bulk) 

7 M \R?p N ) ' 7 ~ M \R 3 pn) ' [ ' 

where S's are dimensionless coefficients depending only on the mode eigenfunctions: 

£(^r) = 1 ( R r 2 drp 2 F(r), vj(^) = I r 2 dr P 2 G{rf . (8.15) 
2 Jo 2 J Q 

Here all quantities in the integrands are in units such that M = R = 1. 

The values of ^ shear "> and Z( bulk ) are given in Table 1 for various modes. Clearly, Z( shear ) increases 
as the order of the mode increases. For higher-order g-modes, our numerical results yield an approximate 
scaling relation, $]^ ft ' eor ) oc u~ 2 . 

8.2 Viscous Dissipation in Coalescing Binaries 

With the results of §8.1, we can now write down the expressions for the viscous energy dissipation rate 
in the coalescing binary neutron stars. Using equations (2.6), (2.16) and (8.3)-(8.5), we have 

E msc = ^(M'Wi m Q n i) 2 \b a - imm a \ 2 2 la = ^{M'Wi m Q n i) 2 \c a - is m uj a c a \ 2 2 la . (8.16) 

a a 

It is clear that the contribution from each mode in this expression is simply twice of the tidal kinetic energy 
(see eq. [6.6]) of the mode divided by energy damping time of that mode, l/(2"f a ). 

Consider first the contribution from the f-mode. In this case, b a is given by equation (3.6). Since 
\b a \ « fi|&a|, we have, after adding contributions from both m = 2 and m = — 2 terms, 

S OTSC , ~2(M'TU ;±2 Qo) 2 |2^ 2 70 , (8.17) 

where we have assumed uo » ft and have used equation (2.4). Here and below, 7 refers to m = ±2 modes, 
thus we have, e.g., j( shear ) = 4! 7 ^ e ar) = Al(j lN R/M){M/p N R 3 ) 2 T l ( - shear '> (see eqs. [8.10], [8.14]). Clearly, 
equation (8.17) agrees with our simple dimensional analysis in (8.2). 
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Now consider the dissipation associated with a specific g-modc. Prior to the resonance, one can still 
use equation (3.6) in (8.16) to obtain an expression similar to (8.17). After the resonance, the tidal energy 
stored in the mode stays nearly constant. The viscous dissipation rate is then given by 

- 2 x m y-r w (8 - l8) 



R yv 100Hz / V - 0003 / \i + q 

where we have used equation (6.11) for AE a = AE n i, and we have included both m = 2 and m = — 2 
contributions. 

The relative importance of the non-resonant viscous dissipation (8.17) and the resonant dissipation 
(8.18) can be compared by calculating the total viscous dissipation accumulated during the inspiral from 
D = oo to the orbital radius before merging D = D m : 



E 



dtE visc — — J -jjt£>E visc , (8.19) 

where tn is given by equation (2.13). For the f-mode non-resonant dissipation, we have 

37r / GM 2 GM \ / i? \ 5 / Rc 2 \ 4 Ql 

m8Cfi " so {— 2l °—) [dZJ [gm) ^ 

,0,3 7? «(^! 270 f), 

where in the second equality we have assumed D m /R = 3 and have used uj 2 ~ 1.5, Qo — 0.56 (see Table 
1). For the g-mode, we can neglect its contribution prior to the resonance, i.e., we assume E visc ^ a ~ for 
D > D a . The total accumulated viscous dissipation due to the resonant g-mode oscillation is then 

.2 / p„2 \ 13/2 / \ 4/3/^7,^2 



5tt 2 _ 7/3 2 (Rx?Y' ( 2 y (GM 2 GM 

E visc , a ~ u82 U a Qniy GM J [l + qj y R c 3 



fm V 7/ V Qm V d3 „-i6/3 ( 2 \ 4/3 {GM 2 GM\ (8 ' 21) 

Because \Q n i\ decreases rapidly as uj n i decreases, only the resonance with the lowest order g-mode is im- 
portant. In obtaining equations (8.20) and (8.21), we have assumed that 7 Q remains constant during the 
inspiral. In reality, the viscosities depend on temperature, which increases as the neutron star is heated up. 
Nevertheless, since ^ hear ) > ^( shear ) ^ e q ua tions (8.20) and (8.21) already indicated that g-mode resonances 
can be important in heating the neutron star compared to the standard f-mode non-resonant heating. 

8.3 Neutron Star Viscosity 

The expressions derived in §8.1 and §8.2 are quite general, and can be used with different viscosities. 
To actually determine the amount of heating of the coalescing binary neutron stars before merging, we need 
to know the viscosity of the nuclear matter. Calculations by Flowers & Itoh (1976, 1979) yield an estimate 
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for the microscopic shear viscosity in liquid core of neutron star. Approximate fitting formulae have been 
given by Cutler and Lindblom (1987). For normal liquid core, the viscosity is dominated by neutron-neutron 
scattering, with 

/ \ 9 ^ 4 

?7 ~ 1.1 x 10 19 T 8 ~ 2 ( — J gem^s -1 , (8.22) 
\PnJ 

where pn = 2.8 x 10 14 g/cm 3 , and T 8 = T/(10 8 K). When the temperature falls below a critical temperature 
<~ 10 9 K, both neutron and proton are likely undergo a phase transition to become a supcrfluid (see, e.g., 
Wambach, Ainsworth & Pines 1991 for review). In this low temperature regime, the viscosity is dominated 
by electron-electron scattering, and it is given by 

V ~ 4.7 x 10 19 T 8 - 2 (^-] gcm-V 1 . (8.23) 
\PnJ 

Since coalescing binary neutron stars are extremely old, their internal heat must have radiated away. All 
cooling calculations indicate that neutron stars older than 10 8 years have core temperature less than ~ 10 6 
K (and surface temperature less than 10 5 K) (e.g., Tsuruta 1992). Therefore, without tidal heating, the 
cores of the neutron stars in the coalescing binary systems are likely to be in the superfluid state. 

The bulk viscosity of neutron star arises from the delay in achieving beta equilibrium as the density is 
changed. Calculation of Sawyer (1989) yields 

C^4.8xl0 18 T 8 6 ^ 2 ( T ^ T )" 2 gcm-V 1 , (8.24) 

where u> is the angular frequency of the perturbation 9 . Because of the strong temperature dependence, bulk 
viscosity is important only when T >> 10 8 K. We shall neglect it in our following discussions. 

Of course, one should be cautious in using these microscopic viscosities. The stable stratification indi- 
cates that turbulent viscosity can not be present in a cold neutron star. Also, tidal distortion is not likely to 
induce turbulence and convection in the binary stars, at least in the linear regime (Seguin 1976). However, 
it is not completely clear whether any other "anomalous" viscosities, such as the mutual friction between 
the electrons and the superfluid neutron vortices (e.g., Sauls 1989, Mendell 1991), can be present in a cold 
neutron star. 

8.4 Heating of the Neutron Stars 

With the microscopic viscosity given by equation (8.23), the mode damping constant is (for m = ±2): 

7«=4x 10- 5 T 8 - 2 £ Q Mi. 4 i?ro 5 s" 1 , (8.25) 
9 This result is based on modified Urea process. When the proton fraction is sufficiently large to allow 
for direct Urea process (Lattimer et al. 1991), the value of £ is larger than equation (8.24) by a factor of 
(p, n /kT) 2 ~5x 10 7 T 8 ~ 2 (where fi n is the Fermi energy of the neutron). Also, equation (8.24) is valid only 
when the the deviation from chemical equilibrium 5/j, = p, n — p, p — p e is much smaller than kT. More general 
expression can be found in Haensel (1992). 
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where we have suppressed the superscript "shear" . 

Below we shall focus on neutron stars with M = 1.4M Q and R = 10 km. Equations (8.17) and (8.25) 
then give the non-resonant heating rate 

T 8 - 2 ergs" 1 , (8.26) 

where we have used S ~ 3 (see Table 1). Similarly, equations (8.18) and (8.25) yield the resonant heating 
rate 

= 13 * 1042 (wk)" 3 (iT&)Vs. e,g,-. ( 8 .2T) 
The heat content of the neutron star is mainly given by the thermal energy of the nonrelativistic 
degenerate free neutron gas in the core: 

U ~ ^-Nk B T^L ~ 4.5 x 10 45 T 2 erg, (8.28) 

where N is the total number of nucleons and ks is the Boltzmann constant 10 . The thermal evolution 
equation of the neutron star during the orbital decay can be simply written as 

— = Evisc + E coo i. (8.29) 
at 

The cooling term E coo i due to neutrino emission and surface photon emission can be shown to be small at 
this low temperature and will be neglected. If we include only the f-mode non-resonant heating E V i SC fi, and 
use equation (8.28), then integrating equation (8.29) yields: 

/QD\ 5/4 

T 8 ~ 0.36 ( . (8.30) 

On the other hand, if we include only the resonant dissipation, we obtain, for D « D n i, 

T 8 ~ 0.42 (-^O 2 (BA. Y'* ( — . (8.3i) 

8 Vioo RzJ Vo.ooo3y V 20 J 17 

Thus the resonant excitation of the g-mode can be as important as the standard f-mode tidal dissipation 
in heating the neutron star. Moreover, this resonant excitation leads to significant heating at larger orbital 
separation, where the f-mode tidal heating is small. Our value of the temperature is much smaller than that 
estimated by Meszaros & Rees (1992). To obtain their value of ~ 10 11 K, we would require a viscosity that 
is larger than the microscopic value by a factor of order 10 10 — 10 12 , corresponding to a kinematic viscosity 



10 When both neutron and proton become superfluid, their heat capacity is greatly reduced by a factor of 
order e A / fcT , where A ~ 0.1 — 1 McV is the superfluid gap energy. In this regime, the main thermal content 
of the neutron star is due to that of the relativistic degenerate electrons, and U is smaller by a factor 2x. 
Using equation (4.25) as a simple estimate of x, we have f/~ 1.2 x 10 44 T| (erg). If we use this expression 
instead of eq. (8.28), the resulting temperature is higher than eqs. (8.30)-(8.31) by a factor of 2.4. 
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comparable to the maximum possible value in a neutron star v max ~ (GMi?) 1 / 2 <~ 10 16 cm 2 /s. We consider 
this rather unlikely. We therefore expect that before the binary neutron stars merge, they remain relatively 
cold and degenerate. Tidal dissipation heats up the neutron stars from less than 10 6 K to ~ 10 8 K before 
the stars come into contact. 

9. CONCLUSIONS 

We have studied in this paper various consequences of tidal excitations of oscillations in coalescing 
binary neutron stars. The energy transfer and angular momentum transfer from the orbit to the star due 
to the resonant excitations of low-frequency g-modes are found to be small because of the weak coupling 
between the g-mode and the tidal potential. The induced orbital phase errors are negligible for the detection 
of gravitational wave from such binary systems by LIGO and VIRGO. Therefore, only the f-modc, quasi- 
cquilibrium tide needs to be considered. Binary models based on ellipsoidal figures (e.g., Lai, Rasio & Shapiro 
1994a, 1994b) are thus expected to give a good description of the effects of tidal interaction on the orbits of 
coalescing neutron star binaries. However, we find that the resonant excitation of the lowest-order g-mode 
can play a significant role in the tidal heating of the neutron stars prior to merger. 

Our calculations of the core g-mode oscillations of neutron stars indicate that the g-mode frequencies 
depend not only on the pressure-density relation, but also on the symmetry properties of the nuclear matter. 
Since the resonant amplitude of the tidal distortion associated with the lowest-order g-mode can be as large 
as ~ 0.5%, if some electromagnetic signature of this induced oscillation can be detected, as suggested by 
RG2, then we may be able to constrain the properties of nuclear matter in the neutron star. We caution 
that the actual values of the g-mode frequencies obtained in this paper are by no means conclusive. In our 
calculations, we have made several assumptions, e.g., we have used Newtonian theory, we have neglected the 
crustal density discontinuities and shear modulus, the thermal effects, as well as possible superfluid effects. 
Also, our calculations are based on a single type of nuclear equations of state (WFF), since these are the only 
EOS available to us from which we can extract enough information to calculate the g-modes self-consistently 
(see §4) . It is desirable that future microscopic calculations of high density nuclear matter not only provide 
the pressure-density relation, but also provide sufficient details which would allow one to obtain the sound 
speed c s and the small difference c 2 — c 2 , so that g-mode properties based on a wider variety of EOS can be 
calculated. 

Finally, we note that although the focus of this paper has been neutron star binaries, the expressions we 
have established in §2-3 and §6-8 are generally valid and can be applied to some other binary systems as well, 
as long as the tidal distortion remains small and the linear approximation is valid. For example, coalescing 
white dwarf binaries have been considered as important sources of low-frequency gravitational waves that 
are potentially detectable by future space-based interferometers (Evans, Iben & Smarr 1987). The restoring 
force for the g-modes of white dwarfs is provided by the thermal energy of degenerate electrons (Baglin & 
Heyvaerts 1969, Chanmugam 1972). Since the orbital decay rate is much smaller for a white dwarf binary, 

31 



the effects of g-mode resonances are much more pronounced. Indeed, using the typical value of the ratio 
Rc 2 /(GM) ~ 10 3 — 10 4 for white dwarfs, and using similar dimensionless values for u n i and Q n i (e.g., 
ujni ~ 0.1, Qni ~ 10~ 4 ), we find that the energy transferred to the white dwarf during a resonance can be 
significant compared to the orbital energy (see eq. [6.13]). Thus g-mode resonances are important for the 
orbital evolution of coalescing binary white dwarfs. 

ACKNOWLEDGEMENTS 

It is a pleasure to thank Curt Cutler, Sam Finn, Edwin Salpeter, Stuart Shapiro, Saul Teukolsky, and 
especially Phil Nicholson and Ira Wasserman for useful discussions on the subjects related to this paper. 
I also thank Andreas Reisenegger, Stuart Shapiro and Saul Teukolsky for reading an early version of this 
paper and helpful suggestions. I am particularly grateful to Andreas Reisenegger for sending his paper to me 
prior to publication and making many critical comments on the earlier version of this paper. Communication 
with Bob Wiringa concerning nuclear equation of state is also greatly appreciated, as well as Greg Cook for 
sharing computer data files of some equation of states for neutron stars. This work has been supported in 
part by NSF Grant AST 91-19475 and NASA Grant NAGW-2364 to Cornell University. 

REFERENCES 

Abramovici, A., et al. 1992, Science, 256, 325 

Arnctt, W.D., & Bowers, R.L. 1977, ApJS, 33, 415 

Baglin, A., & Hcyvaerts, J. 1969, Nature, 222, 1258 

Baym, G, Bethc, H.A., & Pethick, C.J. 1971, Nucl. Phys. A175, 225 

Baym, G, Pethick, C.J., & Sutherland, P. 1971, ApJ, 170, 299 

Bildsten, L., & Cutler, C. 1992, ApJ, 400, 175 

Bradaschia, C, et al. 1990, Nucl. Instrum. & Methods, A289, 518 

Chanmugam, G 1972, Nature, 236, 83 

Cowling, T.G 1941, MNRAS, 101, 367 

Cox, J. P. 1980, Theory of Stellar Pulsation (Princeton University Press) 

Cutler, C, & Lindblom, L. 1987, ApJ, 314, 234 

Cutler, C. et al. 1993, Phys. Rev. Lett., 70, 2984 

Evans, C. R., Iben, I., & Smarr, L. 1987, ApJ, 323, 129 

Finn, L.S. 1987, MNRAS, 227, 265 

Flowers, E., & Itoh, N. 1976, ApJ, 206, 218 

. 1979, ApJ, 230, 847 

Goldrcich, P., & Nicholson, P. 1989, 342, 1079 
Haenscl, P. 1992, A&A, 262, 131 

32 



Higgins, T.P., & Kopal, Z. 1968, Astro. Space Science, 2, 352 

Jackson, J.D. 1975, Classical Electrodynamics, 2nd Ed. (John Wiley & Sons: New York), p. 99. 
Kochanek, C. S. 1992, ApJ, 398, 234 

Kopal, Z. 1978, Dynamics of Close Binary Systems (D. Reidel Publishing Company: Dordrecht-Holland), 
p. 126. 

Lagaris, I.E., & Pandharipande 1981, Nucl. Phys., A369, 470 

Lai, D. 1994, Ph.D. thesis, Cornell University 

Lai, D., Rasio, F.A., & Shapiro, S.L. 1993, ApJ, 406, L63 

. 1994a, ApJ, 420, 811 

. 1994b, ApJ, submitted 

Landau, L.D., & Lifshitz, E.M. 1987, Fluid Mechanics, 2nd Ed. (Oxford: Pergamon Press) 
Lattimcr, J.M., Pcthick, C.J., Prakash, M., & Hacnsel, P. 1991, Phys. Rev. Lett., 66, 2701 
Ledoux, P. 1974, in Stellar Instability and Evolution (IAU Symposium No. 59), ed. P. Ledoux, A. Noels 

and A.W. Rodgers (D. Reidel Publishing Company: Dordrecht-Holland) 
Lee, H.M., & Ostriker, J.P. 1986, ApJ, 310, 176 

McDermott, P.N., Van Horn, H.M., & Scholl, J.F. 1983, ApJ, 268, 837 

McDermott, P.N., Van Horn, H.M., & Hansen, C.J. 1988, ApJ, 325, 725 

Mcndcll, G. 1991, ApJ, 380, 530 

Meszaros, P., & Rees, M. J. 1992, ApJ, 397, 570 

Narayan, R., Piran, T., & Shemi, A. 1991, ApJ, 379, L17 

Nicholson, P. 1978, Ph.D. thesis, Caltcch 

Phinney, E. S. 1991, ApJ, 380, L17 

Prakash, M., Ainsworth, T.L., & Lattimcr, J.M. 1988, Phys. Rev. Lett., 61, 2518 

Press, W. H., Tcukolsky, S. A., Vetterling, W. T, & Flanncry, B. P. 1992, Numerical Recipes: The Art of 

Scientific Computing, 2nd Ed. (Cambridge: Cambridge Univ. Press) 
Press, W.H., & Teukolsky, S.A. 1977, ApJ, 213, 183 
Reisenegger, A. 1994, Preprint 

Reisenegger, A., & Goldrcich, P. 1992, ApJ, 395, 240 (RG1) 
. 1994, ApJ, in press (RG2) 

Sauls, J. A. 1988, in Timing Neutron Stars, eds. H. Ogelman & E.P.J, van den Heuvel (Kluwer: Dordrecht), 
457 

Sawyer, R. F. 1989, Phys. Rev. D, 39, 3804 
Schutz, B. F. 1986, Nature, 323, 310 
Seguin, F. H. 1976, ApJ, 207, 848 

Shapiro, S. L., & Teukolsky, S. A. 1983, Black Holes, White Dwarfs, and Neutron Stars (New York: Wiley) 
Shibata, M. 1993, Osaka University preprint 

33 



Thornc, K. S. 1987, in 300 Years of Gravitation, eds. S. W. Hawking & W. Israel (Cambridge: Cambridge 
University Press), 330 

Tsuruta, S. 1992, in Physics of Isolated Pulsars, eds. K. Van Riper, R. Epstein & C. Ho (Cambridge 
University Press) 

Wambach, J., Ainsworth, T. L., & Pines 1991, in Neutron Stars: Theory and Observations, eds. J. Ventura 

& D. Pines (Kluwer: Dordrecht) 
Wiringa, R. B., Fiks, V., & Fabrocini, A. 1988, Phys. Rev. C, 38, 1010 
Zahn, J.P. 1970, A&A, 4, 452 
. 1977, A&A, 57, 383 



34 



FIGURE CAPTIONS 



Figure 1. — Dimcnsionlcss tidal amplitude associated with a particular g-mode of coalescing binary neutron 
stars as a function of the binary separation D. The resonance radius is assumed to be D a /R = 8, and 
Rc 2 /(GM) = 5, q — M'/M = 1, where R is the stellar radius, M is the neutron star mass, (a) shows the 
real part of the function b a (t) (light solid line) and |6 Q (t)| (dark solid line) for the m = 2 mode; the dashed 
line corresponds to the quasi-equilibrium solution (eq. [3.2]). (b) shows the real part (solid line) and the 
imaginary part (dotted line) of the function c a (t) for m = 2. The orbital phase function is chosen such that 
$(£> = 0) = 0. 

Figure 2. — Equations of state for four neutron star models considered in this paper: (a) pressure P; (b) 
adiabatic sound speed c s (where c is the speed of light); (c) the proton fraction x = n p /n; (d) the fractional 
difference between c 2 s and <? e . The dotted lines are for model AU, the solid lines for model UU, the short- 
dashed lines for model UT, and the long-dashed lines for model UU2. Note that model UU2 has the same 
pressure and c e as model UU. 

Figure 3. — The frequencies of the first ten quadrupole (I = 2) g- modes based on four different neutron star 
models, n specifies the radial order of the g-mode. The dotted line and round circles are for model AU, the 
solid line and squares for model UU, the short-dashed line and triangles for model UT, and the long-dashed 
line for model UU2. 
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